Bulk superconductivity and disorder in single crystals of LaFePO 
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We have studied the intrinsic normal and superconducting properties of the oxypnictide LaFePO. 
These samples exhibit bulk superconductivity and the evidence suggests that stoichiometric LaFePO 
is indeed superconducting, in contrast to other reports. We find that superconductivity is indepen- 
dent of the interplane residual resistivity po and discuss the implications of this on the nature of the 
superconducting order parameter. Finally we find that, unlike T c , other properties in single-crystal 
LaFePO including the resistivity and magnetoresistance, can be very sensitive to disorder. 

PACS numbers: 74.25.Fy, 74.25. Ha, 74.70.-b, 72.80.Ng 



I. INTRODUCTION 

High-temperature superconductivity in the iron- 
pnictides has recently attracted considerable attention 
because it is the first family of materials since the discov- 
ery of the cuprates to exhibit critical temperatures above 
40K. However, in the initial flurry of activity, single- 
crystal measurements of the lower T c analogue LaFePO 
have been somewhat overlooked. Early measurements of 
polycrystalline samples indicated T c values in the range 
from 4 to 7 Kfii2j2iii£ with F-substitution (i.e. electron 
doping) 2 - or Ca-substitution (hole doping) 3 modestly in- 
creasing T c values. However, more recent measurements 
have queried whether the stoichiometric compound is ac- 
tually superconducting, suggesting rather that oxygen 
deficiency plays a significant role£ Initial single-crystal 
work revealed superconducting transitions in the resis- 
tivity and susceptibility, but not in the heat capacity, 
leading the authors to speculate that the superconduc- 
tivity might be associated with an oxygen deficient sur- 
face layer, while the bulk remains metallic to the lowest 
temperatures^ Given the close structural and electronic 
similarity to the parent compounds of the higher T c oxy- 
arsenides, it is particularly important that a clear picture 
of superconductivity in this lower- T c analog is elucidated. 

LaFePO is interesting for reasons other than this con- 
troversy. Because of its lower T c and H C 2 , quantum oscil- 
lations can be measured in relatively moderate fields^ In 
addition, despite the appearance of a nearly-nested Fermi 
surface, LaFePO does not appear to exhibit any mag- 
netism and so the properties of the metallic ground state 
of the pnictides can be examined in detail. In contrast 
LaFeAsO has a spin-density-wave ground state at low 
temperature and only exhibits superconductivity when 
doped, which naturally introduces disorder, inhibiting 
such detailed probes of the Fermiology. 

In this paper we present results of heat capacity, mag- 
netization, resistivity and magnetoresistance measure- 
ments of single-crystal samples of LaFePO. The evidence 
suggests that these samples are stoichiometric, and we re- 
port three major results. Firstly, we find that supercon- 
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FIG. 1: (Color online)(a) Specific heat (left axis) of a mosaic 
of single-crystals and resistivity (right axis) of a single-crystal. 
The heat capacity anomaly corresponds to the disappearance 
of resistance, (b) The susceptibility of powdered LaFePO 
crystals measured in an applied field of 50 Oe following zero 
field cooling (ZFC) and field cooling (FC) cycles, (c) Pow- 
der x-ray diffraction pattern (black) of crystals taken from 
the same batch as in (b). All peaks can be indexed to the 
calculated pattern (red). 



ductivity is present in the bulk of the material, contrary 
to the findings of some other authorsi&i Second, we find 
the superconducting gap is uncorrelated with the inter- 
plane residual resistivity po- If this quantity is represen- 
tative of disorder, the lack of T c suppression is consistent 
with s-wave superconductivity^ Finally, in contrast to 
the superconductivity itself, properties near T c can vary 
greatly with disorder. We illustrate this with measure- 
ments of the magnetoresistance. 
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II. EXPERIMENT 

Single-crystals of LaFePO were grown from a molten 
Sn flux, using conditions modified from the original work 
by Zimmer and coworkers^ Elemental and oxide pre- 
cursors were placed in alumina crucibles and sealed in 
quartz tubes under a small partial pressure of argon. The 
growths were ramped up in temperature to 1190°C fol- 
lowed by a slow cool. We used either La2C>3 or Fe2C>3 
to introduce oxygen into the melt. When La2C>3 was 
used, the highest yield was obtained for a melt contain- 
ing La:La2C>3:Fe:P:Sn in the molar ratio 1:1:3:3:57, while 
with Fe2 03 the best yield was obtained from a combi- 
nation La:Fe2C>3:P:Sn in the molar ratio 3:1:2:24. The 
former method yielded our largest crystals (up to 0.7mm 
on a side) and with the highest RRR (residual resistivity 
ratio, p(300K)/p(0)), but the yield was very sensitive to 
the packing in the crucible, the temperature profile of the 
initial warm up, and the proportion of the other elements. 
The temperature was raised from room temperature to 
1190°C in 6 hours, held for a further 6 hours and then 
lowered to 900°C at a rate of 4.8 °C/hour. Using Fe 2 3 
yielded crystals of the right phase more consistently, al- 
though smaller (typically 0.1-0. 4mm) and with a lower 
RRR. In this case, optimal results were obtained by heat- 
ing over the same amount of time, but the mixture was 
held at 1190°C for 18hrs, before cooling at 10 C/hour to 
650°C. The remaining liquid was decanted at the base 
temperature using a centrifuge, and the crystals removed 
from the crucible. Excess Sn on the surface of the crys- 
tals was removed by etching in dilute HC1 followed by 
an immediate rinse in methanol. Both methods resulted 
in multiphase growths, the principle second phase being 
LaFe2P2, which has the ThCr2Si2 structure. Care must 
be taken to distinguish these two phases. 

Crystals grown using La203 precursor as described 
above had RRR values up to 85, T c ~ 5.9K and an in- 
terplane residual resistivity of po = 0.04mf2cm. Crys- 
tals grown by this technique have been used for re- 
cent ARPES 13 measurements while crystals grown us- 
ing Fe203 were used to measure de Haas- van Alphen 
(dHvA) magnetization oscillations^ The second method, 
using Fe203, with a melt composition much closer to that 
used in Reference 0, yielded crystals with a RRR^ 25, 
T c = 6.7 K and an interplane residual resistivity of 
Po =0.18mf2cm. 

While optimizing the synthesis conditions, batches of 
crystals extracted from individual growths were checked 
by powder x-ray diffraction using an SSI XPERT diffrac- 
tometer. Individual crystals were screened by measure- 
ment of the c-axis lattice parameter. A well-formed single 
crystal grown with La203 precursor and the melt compo- 
sition and temperature profile that produced the highest 
quality crystals (as measured by the residual resistivity) 
was selected for X-ray study. A tin deposit on the crys- 
tal was removed by rinsing in a drop of 12 N HC1. A 
multiply-redundant MoK Q data set of 1643 reflections 
(2&max = 63.9°) was measured on a Bruker Apex II 



diffractometer. The crystal temperature was 90±1 K. 
Data were corrected for absorption vith the face-indexed 
procedure of SADABSl . The structure refined with 
SHELXiS to a standard R index of 0.0084, and a weighted 
R index of 0.02 for all 156 unique reflections. The oxygen 
occupancy is 100% within the statistical uncertainty of 
1.8%. The crystals are tetragonal, space group Ynmm, 
with cell dimensions a = 3.941 ±0.002, c = 8.507 ±0.005. 
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FIG. 2: (Color online) The specific heat C p /T as a function of 
T 2 for a mosaic of LaFePO crystals taken in zero field (solid 
circles) and in an applied field of 1 T (open circles) with field 
oriented parallel to the c-axis. The dotted line indicates the 
linear fit as described in main text which extrapolates to a 
value of 7 ~ 7.0 ± 0.2 mj/molit" 2 . 

Inter-plane (c-axis, p c ) and in-plane (a6-plane, p a b) re- 
sistivity measurements were made with a standard four- 
probe configuration and metallic contact was made using 
Epotek H20E paste. The paste was annealed in air at 
200°C for 20 minutes. Inter-plane measurements were 
made on crystals as small as 0.15x0. 1x0. 02mm 3 and in- 
plane measurements were made on larger crystals. How- 
ever, the larger crystals tended to suffer from inhomo- 
geneity (evidenced by a broad superconducting transi- 
tion), and/or flux inclusions, and the c-axis data were 
generally better. Resistivity and heat capacity mea- 
surements were performed on a Quantum Design PPMS 
(14T). Due to the small size of the crystals, the heat ca- 
pacity was measured for a mosaic of around 100 crystals 
grown from a melt containing La203, with a total mass 
of 800±50/j,<7, each oriented such that the a6-plane was 
parallel to the platform. Crystals from a separate batch 
were powdered, half of which was checked for phase pu- 
rity by powder x-ray diffraction ( Figure[T](c)). Suscepti- 
bility measurements were made for the other half of this 
powder sample using a Quantum Design MPMS 5 Super- 
conducting Quantum Interference Device magnetometer 
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FIG. 3: Heat capacity data on a single-crystal (~5/xg) mea- 
sured using the AC technique giving the relative change of the 
the heat capacity C from the normal state Cn (in arbitrary 
units), where AC = C — C jy . In (a) we show the field depen- 
dence of the heat capacity anomaly, approximately agreeing 
with the suppression of T c extracted from resistivity (see Fig- 
ure [TO]). In (b) we show data fitted with weak-coupling BCS 
theory assuming an isotropic s-wave order parameter. The 
data agrees very well to the theory. 



FIG. 4: (Color online) (a) Interlayer (c-axis) resistivity for 
LaFePO crystals grown using Fe2 03 (blue and green curves) 
and La2 03 (all other curves). The room temperature val- 
ues are in approximate agreement for either method but the 
residual resistance po shows significant variability (inset) . (b) 
In-plane (afe-plane) resistivity. The room temperature value 
shows more variability than the interplane resistivity, primar- 
ily due to difficulties in finding homogeneous samples large 
enough to mount for in-plane measurements. This was often 
manifested in a broader superconducting transition (inset) . 



in a range of applied fields. 

Measurements of heat capacity were also made on a 
single crystal using an AC method.— The sample with 
dimensions 140x220x30p:m 3 (with mass ~ 5/ig) we at- 
tached with GE varnish to a 12/im diameter Chromel- 
Constantan thermocouple, and heated with light at a 
frequency of 16Hz. The size of the temperature oscilla- 
tions is inversely proportional to the heat capacity of the 
sample plus the addenda (varnish plus thermocouple). 
Although with this method it is possible to measure a 
very small single-crystal it is not possible to determine 
the absolute value of the specific heat. The heat capacity 
was measured as a function of field from -0.4T to IT with 



B || c, in a 14T superconducting magnet. To account for 
the residual field in our magnet the actual field was de- 
termined from the field which produced the maximum 
field in a T c versus H plot (the offset was 0.065 T). 



III. RESULTS 

Heat capacity data for the mosaic of single-crystal sam- 
ples described above are shown in Figure [Ua) together 
with representative resistivity data for a single-crystal 
from a separate batch. The total heat capacity of this 
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mosaic is extremely small due to the small mass, making 
measurements of the relaxation time particularly chal- 
lenging. Despite these difficulties, a clear superconduct- 
ing anomaly is evident at 5.9 K, corresponding to the 
temperature at which the resistivity disappears. The 
superconducting state can be suppressed by an applied 
field greater than /xo-ff C 2, leaving only the normal state 
contribution to the heat capacity (Figure \2§ . Data be- 
tween approximately 5 and 10 K can be fitted using 
C v /T = j3T 2 + 7 where (3 and 7 are the lattice and elec- 
tronic heat capacity coefficients respectively. This yields 
an estimate of 7 = 7.0 ± 0.2 mJ/moLKT 2 where one mole 
refers to the formula unit LaFePO. This value is similar 
to that obtained from values extracted by other authors 
for polycrystalline^i^ and single-crystal^ samples, though 
care should be taken on the definition of a mole (the value 
is doubled if one mole refers to the unit cell). Data on the 
low-temperature side of the superconducting transition 
are less reliable due to the small value of the heat capac- 
ity and the experimental difficulties measuring the relax- 
ation time as described above. Due to these difficulties, 
the data does not appear to conserve entropy. Combined 
with the rounding of the superconducting anomaly seen 
in the data, it is difficult to obtain an accurate estimate 
of the heat capacity jump, but we can obtain a conserva- 
tive lower bound by comparing the peak of the transition 
with the normal state extrapolation (vertical arrow in 
Figure [T]). This results in an estimate of the normalized 
heat capacity jump AC/jT c w 0.6 ±0.2, which is greater 
than Kamihara et aL's estimate based on measurements 
of polycrystalline samples^ and is less than (about 40%) 
the BCS result (AC/-/T c ) BC s = 1-43. The actual value 
is likely somewhat larger than this lower bound. Suscep- 
tibility measurements for a sample consisting of several 
powdered single-crystals are shown in Figure [T](b) for an 
applied field of 50 Oe following both zero field cooling 
(ZFC) and field cooling (FC) cycles. Perfect diamag- 
netism (AttM/H = — 1, shown on the right axis in Figure 
lb) corresponds to a value of -1.32x 10~ 2 emu/ g. Assum- 
ing that the particle size is larger than the penetration 
depth (~ 0.2/zro±l), the FC value ( 0.32emu/g) pro- 
vides the better estimate of the superconducting volume 
fraction, giving a lower bound of approximately 24%. 
This likely substantially underestimates the actual value 
due to flux pinning effects. 

Heat capacity data for a single crystal, obtained using 
the AC technique, are shown in Figure [3l For a field of 
> 0.47T the superconducting anomaly was suppressed to 
below T=2.3 K so we used this data as an estimate of the 
normal state heat capacity Cn (including addenda), and 
subtracted this from the data in zero field. The supercon- 
ducting anomaly is now clearly visible, with T c (mid point 
of the increase in C of T=5.6 K and width (10-90%) 
of 0.4K. In Figure [3] (a), we show the field dependence 
of the anomaly in C. The decrease in T c corresponds 
to dT c /dH = 7.2 ± 0.1K/T. Also shown Figure [3] (b) is 
the temperature dependence calculated from the weak- 
coupling BCS theory assuming an isotropic s-wave order 
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FIG. 5: Plot showing the dependence of the critical temper- 
ature T c on the residual resistivity po extracted from the low 
temperature normal state of LaFePO. The higher residual re- 
sistivities correspond to a ratio of the coherence length to 
mean free path £ a b/l ~ 0.17 — 0.25, and so a partial, sys- 
tematic suppression in T c is expected for an anisotropic order 
parameter. 
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FIG. 6: (Color online) (a) c-axis resistivity data from Figure 
3(a) replotted on a logarithmic scale. Plotted this way, all 
data fall on the same curve. The black line shows the value 
A ~ 3.2 x 10~ 5 nrs7cm/i(" 2 , while the grey lines illustrate our 
error bars in A. (b) The coefficient A extracted from fits of 
the resistivity data to p — po + AT 2 , which is independent of 
disorder. 



parameter. Apart from the data very close to T c or for 
T<2.8K the theory fits the data very well. The departure 
below ~2.8K probably arises from the field dependence 
of the thermocouple thermopower and field dependence 
of the addenda which we have not attempted to correct 
for. 

In the strong coupling Eliashberg theory the slope near 
T c normalised to the jump at T c R = 3^§y |t c is indica- 
tive of the coupling strength^ The fact that our data fit 
well the weak coupling BCS form shows that LaFePO is 
in the weak coupling limit (note however this does not 
prove the order parameter is necessarily isotropic). 

In FigureUfa) we show the inter-plane (c-axis) resistiv- 
ity p c for five crystals grown from slightly different melt 
conditions. Absolute values of the room temperature re- 
sistivity are in excellent agreement, and the principal dif- 
ference between the samples is the residual resistivity po . 
A small upturn is seen in the resistivity of some but not 
all samples, just above T c . T c values are also in excel- 
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FIG. 7: (a) The magnitude of the interplane resistivity up- 
turn Ap c — p p — po normalized to the residual resistivity po, 
where p p is the peak resistivity of the upturn. There does not 
appear to be a correlation with T c . (b) The magnitude of the 
interplane resistivity upturn Ap c as a function of the inter- 
plane residual resistivity po, which also shows no correlation. 



lent agreement, and show no correlation with the resid- 
ual resistivity (see Figure [5]). Similar data for in-plane 
resistivity measurements p a b are shown in Figure HJb) . 
In contrast to the c-axis data, in-plane measurements 
tend to exhibit both broader transitions and also a wider 
variation in the absolute value of the resistivity. As de- 
scribed above, this difference can be ascribed to the use of 
larger crystals for the in-plane measurements, which are 
consequently more susceptible to defects, flux inclusions 
and, due to the malleability of the crystals, macroscopic 
distortion. The anisotropy of the interplane to in-plane 
resistivity p c / p a b is in the range from 13 to 17, in approx- 
imate agreement with band structure calculationsj 18 ' 19 i 20 

Neglecting the low-temperature upturn seen in some 
samples, resistivity data below 25K can be well-fit by 
the standard Fermi liquid behavior, p = po + AT 2 , where 
Po is the residual resistivity. 21 To confirm the reliability 
of this fitting, log(p — po) vs log(T) is shown in Figure 
[(Ja) for the c-axis data, which reveals that each data set 
not only has a gradient of 2 ± 0.3, but also falls on the 
same value of A ~ 3.2 ± 0.5 x 10~ 5 mf2cm/if 2 (Figure 
[SJb)). We note in passing, that when A is combined 
with our value of 7 we find a Kadowaki- Woods ratio of 
A/7 2 = 65 ± 15 x 10" 5 ^cm(mol K/mJ) 2 3 

As mentioned above, a small resistivity upturn is ob- 
served in some, but not all, samples just above T c . This 
effect is seen for both current orientations (i.e. for p a b 
and p c ). Significantly, the size of the upturn is not cor- 
related with either T c (Figure E(a)), or with the abso- 
lute value of the residual resistivity po extrapolated from 
higher temperatures (Figure E[b)). We return to the field 
dependence and origin of this feature below. 

Representative data showing the variation in the c- 
axis resistivity as a function of applied field at various 
temperatures are shown in Figure [5] for fields applied 
(a) parallel and (b) perpendicular to the c-axis. From 



data like these we are able to extract p.oH c2 , choosing 
for simplicity to define p-oH C 2 in terms of the midpoint 
of the resistive transition, and using 10 and 90% of the 
transition to determine the error bars. Plotting p.oH C 2 
extracted in this manner, as shown in Figure [9] and ex- 
trapolating to T = OK, we determine pqH^O) ~ 0.6T 

and pqH^ 2 (0) ~ 3.5T, giving an anisotropy of ~ 6. There 
is some variation in this value, perhaps due to the sen- 
sitivity of _ffj 2 to orientation, but most samples yield an 
anisotropy between 6 and 15, similar to that observed by 
other authors^ 

The resistivity upturn which is observed in several 
samples is suppressed by an applied field, resulting in 
a pronounced "hump" feature in the magnetoresistance 
p(H) for fields just above H C 2 (see Figure [5] ). The 
magnitude of the resistivity upturn in zero field, and 
hence the associated negative magnetoresistance, varies 
between samples grown under different conditions, and 
to a lesser extent even between samples from the same 
batch. This variation is illustrated in Figure [TU] which 
show representative c-axis and in-plane magnetoresistiv- 
ity data for samples with small (panel (a)) and large 
(panel (b)) upturns respectively. The field is oriented 
parallel to the c-axis in both cases to allow direct com- 
parison. For the sample shown in panel (a), the anomaly 
is suppressed by fields as small as 0.3T, whereas for the 
sample shown in panel (b) it survives to fields of more 
than twice this value. In all cases, once the anomaly 
is completely suppressed, a positive magnetoresistance 
appears which is approximately linear with field in the 
high- field limit. 

Significantly, the suppression of the resistivity upturn 
is strongly dependent on the orientation of the applied 
field with respect to the c-axis. For example, for the 
sample shown in Figure [TU (a-c), the resistivity upturn 
is suppressed for fields exceeding a modest 1 T when ap- 
plied parallel to the c-axis, but is only suppressed for 
fields exceeding 12 T when directed perpendicular to the 
c-axis. Although the absolute values of the fields at which 
the resistivity upturn are suppressed vary between crys- 
tals depending on the magnitude of the zero-field upturn, 
this anisotropy is always observed. This anisotropy ap- 
pears to mirror the anisotropy in poH C 2, such that if 
magnetoresistance data for the two orientations are plot- 
ted as a function of field scaled by the appropriate value 
for p.oH C 2 for that orientation, the data lie almost on top 
of each other (Figure [TT1 (d)). It is not clear whether 
this is coincidental, or reflects a deeper inter-relation be- 
tween superconductivity and the origin of the resistivity 
upturn. 



IV. DISCUSSION 

The results summarized in Figure [Tj3] provide the first 
evidence of superconductivity in single-crystals of sto- 
ichiometric LaFePO, manifested in a bulk thermody- 
namic transition in both the magnetization and the heat 
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FIG. 8: (Color online) Interlayer (c-axis) resistivity as a func- 
tion of magnetic field for (a) field parallel and (b) orthogonal 
to the c-axis for temperatures from 1.8 to 6.0 K (see legend 
in upper panel). The anisotropy in H C 2 is approximately 6 in 
this case, though other samples show an anisotropy of up to 
15. 



capacity. In contrast, despite the observation of bulk 
superconductivity in polycrystalline samples by several 
other groups^£ McQueen et al. have recently sug- 
gested that the stoichiometric compound is actually non- 
superconducting, and that the onset of superconductivity 
is perhaps associated with the presence of oxygen de- 
ficiencies, among other possibilities^ However, the ab- 
sence of superconductivity in LaFePO would be surpris- 
ing given that there is nothing anomalous in the density 
of states^ of the stoichiometric compound, no evidence 
of a competing phase transition, and doping with either 
holes or electrons results in superconductivity^ Using 
crystals grown by a similar technique to our own, Hamlin 
et al. observed a substantial diamagnetic susceptibility, 
but no anomaly in the heat capacity, causing them to sug- 
gest the possibility that their crystals have a thin surface 
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FIG. 9: The critical field fioH C 2 for fields oriented parallel 
(solid symbols) and perpendicular (open symbols) to the c- 
axis. Data points were determined from resistivity data as 
described in the main text. Error bars indicate 10 and 90 % 
of the superconducting transition. 



layer of oxygen-deficient superconducting material sur- 
rounding a stoichiometric non-superconducting center.— 
Refinements of single-crystal x-ray diffraction data for 
our crystals indicate full site occupancy, with a standard 
uncertainty in the oxygen content of only 1.8%. This 
is not the best probe of the oxygen content, but several 
other pieces of evidence suggest that oxygen deficiency, 
or at least variation in oxygen deficiency, does not play 
a significant role in these crystals. Specifically, neither 
the T 2 coefficient of the resistivity nor T c itself (Figures 
[6] and [5] respectively), show any significant variation with 
the residual resistivity, in contrast to what would be ex- 
pected if variation in the oxygen content were responsible 
for the different residual resistivity values if this were the 
principle source of disorder in these samples. Further- 
more, these values do not depend on details of the syn- 
thesis conditions, including whether oxygen is introduced 
to the melt using Fe 2 03 or La2C>3, and even the rela- 
tive concentration of these precursors in the melt. These 
observations demonstrate that there is remarkably little 
variation in oxygen deficiency between crystals grown by 
this technique. In addition, the extremely high resid- 
ual resistivity ratios, up to 85, indicate that the material 
is in fact exceptionally well-ordered. Finally, analysis 
of recent dHvA data taken for samples grown by this 
technique indicate that the electron and hole pockets are 
almost balanced, in agreement with the anticipated re- 
sult for stoichiometric material^ with a mismatch that 
would correspond to a doping smaller than the error of 
our x-ray diffraction refinements. 

Experiments involving annealing in an oxygen atmo- 
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sphere may be revealing, and though these are in progress 
one will always have the problem of establishing by how 
much the material has departed from perfect stoichiome- 
try, which is a difficult measurement for such small crys- 
tals. Thus we cannot definitively prove the absence of 
oxygen deficiency from the present measurements, and 
in fact the breadth of the specific heat anomaly in Fig- 
ure [T] may be indicative of slight variation in stoichiom- 
etry even within a single sample. However, given the 
small departure from stoichiomctry that our errors al- 
low, together with the above observations, it seems ex- 
tremely unlikely that the stoichiometric compound is 
non-superconducting. Furthermore, the absence of other 
competing phases suggests that superconductivity is the 
natural low temperature ground state. In this case, the 
absence of a superconducting anomaly in the heat capac- 
ity data of Hamlin et al. likely reflects similar difficulties 
that we experienced in the measurement of such small 
sample masses. It is difficult to speculate on the rea- 
son for the absence of superconductivity in the polycrys- 
talline samples of McQueen et al£ but we note that the 
resistivity of those samples exhibits a substantial low- 
temperature upturn. We return to the possible signifi- 
cance of this below. 

From the range and variety of samples measured in 
this work, we are also able to consider the role of dis- 
order on the superconductivity of LaFePO. Even though 
these materials are anisotropic, the interplane transfer in- 
tegral for the electron pockets is relatively large^ and so 
the interplane residual resistivity should be a good mea- 
sure of the number of defects in a given crystal, though 
this argument is weaker for the hole pockets which have 
little warping in the k z direction^ Assuming that this is 
the case, the plot of T c vs po, as shown in Figure O illus- 
trates that T c is, within our uncertainties, uncorrelated 
with disorder. For other high-temperature superconduc- 
tors the absence of T c suppression has been attributed 
to the smallness of the coherence length 2 ^ t; a b- From the 
present measurements ^ a b can be estimated from the rela- 
tion HoH<k = $o/27r£a 6 , giving £ a b ~ 0.03pm. The mean 
free path can be estimated from a quasi-two-dimensional 
Drude model mi = (e 2 /h) x kpl/c, which yields a dif- 
ferent conductivity for electron and hole pockets. Us- 
ing the results of Reference [H, we estimate the electron 
and hole pockets to have a mean free path I ~ 1700 A 
and I ~ 1200A respectively for the dirtiest samples, ap- 
proximately in agreement with estimates from de Haas 
van Alphen oscillations. 8 For the electron (hole) pockets 
this suggests the ratio £ a b/l has shifted from a value of 
~ 0.04 (0.06) for the cleanest samples to ~ 0.17 (0.25) 
for the dirtiest samples. If the disorder were magnetic, 
a partial, systematic suppression of T c should manifest 
for any pairing symmetry according to the relationship 
first calculated by Abrikosov and Gor'kov i 25 ' 26 Similarly, 
T c for anisotropic ci-wave or p-wave superconductors in 
the weak coupling limit should be suppressed due to the 
inevitable phase mixing of the order parameter, which 
should be noticeable by this value of £/?i 24 i 26 i 27 Con- 



versely, according to Anderson's theorem^ an s-wave 
superconductor should be robust against non-magnetic 
disorder due to the ability of the quasiparticle states to 
form time reversed pairs. The data in Figure suggests 
that the predominant form of disorder is therefore non- 
magnetic. More importantly, if the dominant scattering 
is strong enough to scatter across the Brillouin zone, then 
this plot also illustrates that the order parameter mediat- 
ing superconductivity is likely to be s-wave, an observa- 
tion consistent with recent measurements of penetration 
depthi 7 - on SmFeAsOi_ ;E F. E . However, if the scattering 
is small angle, or only occurs intra-Fermi surface, then 
an s-wave order parameter which reverses sign between 
Fermi surfaces centered at the T and M points may also 
be possible, as proposed by Mazin et al£& On the other 
hand, if the disorder is inhomogeneous on scales of order 
£, then it is possible than an anomalously clean path will 
yield a resistive transition T c which is essentially disor- 
der independent, though the superfluid density would be 
rapidly suppressed. This behavior may be evident in the 
heat capacity, and the breadth of the anomaly shown in 
Figure [1] may be indicative of just such inhomogeneity. 29 
Naturally, the above discussion relies on the interplane 
resistivity being an accurate measure of the in-plane scat- 
tering rate, which may be true for the electron pockets 
due to their significant c-axis dispersion, but less likely 
for the hole pockets. 

We now turn the discussion to the origin of the re- 
sistivity upturn observed in some crystals for temper- 
atures just above T c . Similar upturns are observed in 
several other classes of unconventional superconductors 
including many of the cuprates* 3 - - several organic o 31 i 32 
and other polycrystalline oxy-arsenides^ 3 - However, in all 
these classes of compounds one observes an enhancement 
of the anomaly with an applied magnetic field, whereas in 
the present case we observe a suppression of the anomaly 
with field. In addition, we observe the upturn in some 
but not all of our samples, implying that this behavior is 
intimately linked with the presence of disorder. Finally, 
the anomaly is observed regardless of the current orien- 
tation, whereas in other correlated materials the effect 
is only observed when the current is passed along the 
c-axis. 30 In highly anisotropic materials with Josephson 
coupled planes, one expects dramatic differences in the 
behavior of each current orientation in field. Specifically, 
near T c in-plane currents are dissipated by the motion 
of pancake vortices under the Lorentz force while inter- 
layer currents are dissipated by the Josephson coupling 
between the layers, giving rise to a broad anomaly in 
p c 22- We thus rule out such coupling as the source of the 
anomaly. 

Given the relatively large concentration of Fe in the 
melts used to grow these crystals, it is not unreasonable 
to suggest that the resistivity upturn might be associated 
with Kondo scattering from magnetic impurities. How- 
ever, the strong anisotropy of the upturn's suppression 
with field, together with the fact that samples exhibit- 
ing the upturn do not show a dramatically suppressed 
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FIG. 10: (Color online) Temperature sweeps at constant field 
for field directed parallel to the c-axis illustrating the effect of 
applied fields on the resistivity upturn observed at low tem- 
peratures for some samples. In (a) we show the interplane 
resistivity for a sample with a small upturn and in (b) we 
show the in-plane resistivity for a sample with a pronounced 
up-turn. In both cases the anomaly is rapidly suppressed by 
magnetic field, though clearly the more pronounced the up- 
turn, the higher the field required for the suppression. 



T c convincingly excludes this possibility. Another pos- 
sible scenario for the origin of the resistivity anomaly is 
weak localization due to elastic scattering defects, possi- 
bly also explaining the rapid suppression of the anomaly 
with field. Furthermore, in quasi-two-dimensional sys- 
tems, there can be a significant anisotropy in the field- 
induced suppression of weak localization^ However, the 
upturn does not scale with the residual resistivity po (as 
shown in Figure [7] (b)), indicating that the effect does 
not correlate with the concentration of elastic scattering 
centers. 

An intriguing alternative hypothesis for the origin of 
the resistivity anomaly is scattering from spin fluctua- 
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FIG. 11: (Color online) In-plane magnetoresistance for a 
sample with a pronounced upturn. At low fields the magne- 
toresistance is negative until the anomaly is completely sup- 
pressed, (a) illustrates the suppression when field is applied 
in-plane while (b) shows that the anomaly extends over a 
much wider field range when the field is applied orthogonal 
to the c-axis. (c) Illustrates both these data set on the same 
plot for direct comparison and (d) shows the data at 1.8 K 
scaled by fioH^ for the in-plane data and ^oH^ 2 for the inter- 
plane data. 



tions associated with a competing magnetic phase. Near 
perfect nesting of electron and hole pockets evidenced 
from dHvA measurements^ certainly raises the question 
as to whether the spin density wave (SDW) observed 
in the isoelectronic compound LaFeAsO might be sta- 
bilized in LaFePO. In contrast, recent NMR results for 
Lao.87Cao.i3FePO indicate the presence of ferromagnetic 
fluctuations^ perhaps indicating that q ~ nesting 
between concentric FS sheets dominates for such high 
doping levels. The same study also revealed anomalous 
behavior of 1/T\T below T c for that material, indicat- 
ing unusual spin dynamics in the superconducting state, 
though the origin of this effect is not yet clear. While the 
evidence points towards stoichiometric LaFePO being su- 
perconducting, one can speculate that, for example, mag- 
netic impurities, or even disorder-induced variation in the 
Fe-P bond-length (which is believed to affect the magni- 
tude of the Fe moment 36 ) might stabilize an incipient 
SDW in this material. However, there is no evidence for 
a magnetic or density wave transition in single-crystals 
of LaFePO, and T c is uncorrelated with the magnitude 
of the resistivity upturn, at least for the range of val- 
ues that we have observed. Hence, for such a scenario 
to be applicable, either the superconductivity must be 
insensitive to the presence/absence of magnetic fluctu- 
ations, or there must be some degree of inhomogeneity 
in the crystals, such that parts of the sample that suffer 
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such fluctuations are physically separate from the major- 
ity superconducting phase. Unfortunately although this 
possibility is especially interesting due to the parallels to 
LaFeAsO, we cannot definitively establish whether it is 
appropriate from these measurements. Additional exper- 
iments are underway to determine the effect of stabilizing 
the SDW on the normal and superconducting states in 
this compound. 

We note in closing that whatever the cause of the re- 
sistivity anomaly, it is reminiscent of an upturn seen by 
McQueen et al. for polycrystalline samples which had 
no superconducting transition^ The residual resistivity 
ratio of those polycrystalline samples was nearly two or- 
ders of magnitude lower than that of the single-crystal 
samples described in this report, and the resistivity up- 
turn was significant. This raises the possibility that the 
lack of superconductivity in those samples might not be 
due to differences in the oxygen content, but rather to 
a larger concentration of defects, potentially affecting a 
competing density wave ground state. 

V. CONCLUSIONS 

We have studied transport and thermal properties of 
single-crystals of LaFePO. These samples exhibit super- 



conductivity, and the evidence suggests that stoichiomet- 
ric LaFePO is indeed a bulk superconductor, in contrast 
to other reports. Furthermore we find that the super- 
conducting gap is unrelated to non-magnetic disorder 
for values of £/l ~ 0.2(assuming that disorder is pro- 
portional to the interplane resistivity). This result is 
consistent with an s-wave pairing symmetry, but is not 
definitive. Finally, we have investigated an anomalous 
low-temperature resistivity upturn that appears in some 
of our samples through detailed temperature and field 
studies. The feature is related to disorder, though the 
present measurements are insufficient to definitively es- 
tablish the origin of this effect. 
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